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Abstract. Simulations are performed to investigate turbulent properties of non- 
linearly interacting two-dimensional (2D) magnetic electron drift vortex (MEDV) 
modes in a nonuniform unmagnetized plasma. The relevant nonlinear equations gov- 
erning the dynamics of the MEDV modes are the wave magnetic field and electron 
temperature perturbations in the presence of the equilibrium density and tempera- 
ture gradients. The important nonlinearities come from the advection of the electron 
fluid velocity perturbation and the electron temperature, as well as from the nonlin- 
ear electron Lorentz force. Computer simulations of the governing equations for the 
nonlinear MEDV modes reveal the generation of streamer-like electron flows, such 
that the corresponding gradients in the direction of the inhomogeneities tend to 
flatten out. By contrast, the gradients in an orthogonal direction vary rapidly. Con- 
sequently, the inertial range energy spectrum in decaying MEDV mode turbulence 
exhibits a much steeper anisotropic spectral index. The magnetic structures in the 
MEDV mode turbulence produce nonthermal electron transport in our nonuniform 
plasma. 



1. Introduction 

About a quarter century ago, several authors [1, 2, 3] predicted the existence of 
the magnetic-electron-drift vortex (MEDV) mode in a nonuniform unmagnetized 
plasma because of its possible importance in laser fusion as well as in nonstationary 
processes in laboratory plasmas. The MEDV mode is purely magnetic and has negli- 
gible density and electrostatic potential perturbations in an inhomogeneous plasma 
with fixed ion background. When the equilibrium electron density and electron 
temperature gradients are simultaneously present, the MEDV mode gets unstable 
[2, 3] because of a first order baroclinic force, which can enhance the wave magnetic 
field and electron temperature perturbations if the electron temperature gradient is 
sufficiently large in comparison with the density gradient. Large amplitude MEDV 
modes are subjected to parametric instabilities [4, 5] involving the decay and mod- 
ulational interactions [6, 7, 8]. Due to the parametric processes the magnetic energy 
of the MEDV modes is redistributed among other modes in plasmas. 

In their classic paper, Nycander et al. [9] derived a pair of equations governing the 
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dynamics of nonlinearly interacting MEDV modes. It turns out that the nonlinear 
mode coupling equations contain the intrinsic Jacobean nonlinearities arising from 
the advection of the electron fluid velocity and electron temperature perturbations, 
as well as the nonlinear Lorcntz force. These nonlinearities are responsible for the 
formation of coherent magnetic structures in the form of a dipolar vortex [10, 11] 
and a vortex street [12]. 

Recently, there has been a renewed interest in the study of MEDV mode turbu- 
lence [12, 13]. Specifically, Jucker and Pavlenko [12] and Andrushchenko et al. [13] 
have presented analytical studies of the generation of large-scale magnetic fields 
(magnetic streamers) via the modulational instability of coherent and partially in- 
coherent MEDV modes. The use of a quasi-linear theory also reveals the shearing 
of micro-turbulence by the electron flows and the corresponding electron diffusion 
in the presence of zonal magnetic fields and magnetic streamers [13] . 

In this paper, we present results of computer simulations of 2D nonlinearly in- 
teracting MEDV modes whose dynamics is governed by the nonlinear equations in 
Ref. [10]. Needless to say, we numerically solve the coupled equations for the wave 
magnetic field and the electron temperature perturbation in order to investigate 
the formation of magnetic and temperature fluctuation structures, as well as the 
associated turbulent spectra and the resultant electron transport in the magnetic 
structures. 



2. 2D nonlinear equations 

We consider a nonuniform plasma in the presence of the equilibrium density and 
electron temperature gradients. The latter, which are along the x axis in a Cartesian 
coordinate system, are maintained by the external sources, such as the dc electric 
field . The dynamics of the MEDV mode in our nonuniform plasma is governed 
by the electron momentum equation 

^ -Fu- V ) u-h-^ (E-hiu X B ) -h— V(nT) = 0, (1) 



^dt J m\ c / mn 

the electron energy equation 



^+u-V)r+(7-l)TV-u = 0, (2) 



the Faraday law 
and the Amperes law 



VxE = -lf, (3) 

„ „ AirenVL ,,, 
V X B = , (4) 

c 

where u is the electron fluid velocity, T is the electron temperature, e is the mag- 
nitude of the electron charge, m is the electron mass, c is the speed of light in 
vacuum, 7(= 2/3) is the adiabatic index (the ratio of the specific heat), and E and 
B are the electric and magnetic fields, respectively. The displacement current in (4) 
has been neglected since the MEDV mode phase speed is much smaller than c. The 
assumption of immobile is justified as long as the MEDV mode frequency is larger 
than the ion plasma frequency. 

Taking the curl of (1), using (3) and (4) and letting B = zB, T = To{x) + Ti, 
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where z is the unit vector along the z axis and ri(<C To) is a smah temperature 
fluctuation in the equihbrium value Tq{x), we obtain [10] 

I (B - A^V^B) + —[B,V'B], (5) 

ot ^ 'ay mc 

where /3 = cKn/e, Kn = dlnn{x)/dx, A = c/ujp, and uip = {inne^ /m)^^^ is the 
electron plasma frequency. Wc have denoted = {d^ / dx^) + {d^ / dz^) and [F, G\ = 
z X VF • WF represents the Jacobean nonlinearity. 
The electron energy equation (2) is written as 

§ + ,|? + £^lB,r] = o, (6) 

at oy mc 

where a = \^{eT(,/mc)[{2/i)Kn - Kt], and Kt = |91nTo/aT|. 

In the linear limit, we neglect the nonlinear terms in (5) and (6), and by assuming 
that B and Ti are proportional to exp(— iwf + ik • r), we Fourier transform and 
combine the resultant equations to obtain the MEDV mode frequency 

w = fcj,Va/?/(l + fc2A2), (7) 

where a/3 = V^\^Kn[{'2./i)Kn - Kt] and Vr = (T/mYf^ is the electron thermal 
speed. Furthermore, k"^ = ky + k"^ and ky{kx) is the component of the wave vector 
k along the y{x) axis. The wavelength of the MEDV mode is assumed to be much 
smaller than K""^ and K^"^ . Equation (7) reveals that for fc^A^ <C 1, the phase 
speed uj/ky of the MEDV mode is constant, while for /c^A^ ^ 1, we have lo = 
{ky/k)^Ja(3/X^. The latter indicates frequency condensation at kx ~> ky. 
Equations (5) and (6) admit the energy integral 

E = J[B^ + X^^Bf + {0/a)T^]dxdy. (8) 

The total energy above comprises the magnetic field, magnetic field vorticity and 
temperature energies. The energy integral (8) is used to check the numerical validity 
of our code, which must be preserved under ideal conditions. 



3. Simulation results 

We suitably normalize (e.g. t by uj^^, V by A~^, B by mcojp/e, Ti by Tq ) (5) and 
(6) and investigate numerically the properties of 2D localized magnetic structures, 
as well as examine the associated magnetic turbulence spectrum and nonthermal 
cross-field electron transport caused by the magnetic structures. 

Wc have developed a code to integrate the normalized (as mentioned above) 
equations (5) and (6). The temperature fluctuaion is reinforced due to nonlinear 
couplings between random set of magnetic field fiuctuations and an initial temper- 
ature perturbation, while enhanced temperature fluctuations drive magnetic field 
fluctuations linearly due to the baroclinic force. The nonlinear interactions between 
MEDV magnetic fields keep mainitianing the resulting magnetic structures. Our 
numerical code employs a doubly periodic spectral discretization of magnetic field 
and temperature fiuctuations in terms of its Fourier components, while nonlinear 
interactions are deconvoluted back and forth in real and Fourier spectral spaces. 
The time integration is performed by using the Runge-Kutta 4th order method. 
A fixed time integration step is used. The conservation of energy (8) is used to 
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Figure 1. Formation of streamer- like fiows in nearly steady-state magnetic field fluctua- 
tions. The temperature fluctuation is predominantly nearly isotropic, but some eddies are 
aligned in the direction of the background inhomogeniety. Numerical resolution is 512^, 
box dimension is 2ii x 2ti, a = /3 — 1.0. 

check the numerical accuracy and vaUdity of our numerical code during the nonhn- 
ear evolution of the magnetic field and temperature fluctuations. Varying spatial 
resolution (from f28^ to 512^), time step (10~^, 5 x 10~^, 10~^), constant values of 
KnXV^/c^ = 0.1 and {2/3)KnX — KtX = 0.1 are used to ensure the accuracy and 
consistency of our nonlinear simulation results. We also make sure that the initial 
fluctuations are isotropic and do not influence any anisotropy during the evolution. 
Anisotropy in the evolution can however be expected from & ky = Q mode that is 
excited as a result of the nonlinear interactions between the finite frequency MEDV 
modes. 

The magnetic and temperature fluctuations are initialized by using a uniform 
isotropic random spectral distribution of Fourier modes concentrated in a smaller 
band of wave number (/c < 0.1 kmax)- While spectral amplitude of the fluctuations 
is random for each Fourier coefficient; it follows a k~^ or fc~^ scaling. Note again 
that our final results do not depend on the choice of the initial spectral distribution. 
The spectral distribution set up in this manner initializes random scale turbulent 
fluctuations. Since there is no external driving mechanism considered in our sim- 
ulations, turbulence evolves freely under the influence of self-consistent instability. 
Note however that driven turbulence in the context of the MEDV mode will not 
change inertial range spectrum to be described here. 

The nonlinear evolution of the magnetic and electron temperature fluctuations is 
governed typically by different nonlinearities in the two dynamical equations (5) and 
(6). To understand the characteristics of the nonlinear interactions in our MEDV 
turbulence model, it is interesting to compare it with that of the Hasegawa-Mima- 
Wakatani (HMW) model [14, 15, 16, 17] describing the electrostatic drift waves 
in an inhomogeneous magnetoplasma. For instance, the generalized magnetic field 
vorticity is driven by z x VB ■ W^B nonlinearity. The latter is similar to the ion 
polarization drift nonlinearity z x Vi^ • VV^0, where cj) is the electrostatic poten- 
tial, in a pseudo-3D HMW equation and signifies the Reynolds stress forces that 
play a critical role in the formation of zonal-flows [8]. Analogously, one can expect 
the nonlinear generation of electron flows in our MEDV model. The nonlinearity 
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Figure 2. Formation of relatively smaller scales streamer-like flows in generalized and 
magnetic field vorticity fluctuations corresponding to the magnetic field structures shown 
in Fig (1). As expected, spectra of the latter contain smaller scales with increasingly 
dominated ky ~ modes. 



in the temperature perturbation equation (6) is z x VTi • V-B, which is akin to a 
diamagnetic nonlinear term in the HMW model. The role of this nonlinearity has 
traditionally been identified as a source of suppressing the intensity of the non- 
linear flows in the drift wave turbulence. Nevertheless, the presence of the linear 
inhomogenous background in both the equations can modify the nonlinear mode 
coupling interactions in a subtle manner. Our objective here is to understand the 
latter in the context of the MEDV turbulent processes. The initial isotropic and 
homogeneous spectral distribution in magnetic field and temperature fluctuations, 
as described above, evolve dynamically following the set of equations (5) and (6). 
The spectral magnetic energy from large scale eddies migrates towards the smaller 
eddies following the Richardson's cascade law. In configurational space, this essen- 
tially corresponds to breaking up the larger eddies into the smaller ones. Conse- 
quently, the mode coupling interactions in the Fourier space follow a Kolmogorov- 
type phenomenology in that spectral transfer, which predominantly occurs in the 
local /c-space in the inertial range MEDV mode turbulence. During this process, 
each Fourier modes in the inertial range spectrum obey the vector triad constraints 
imposed by the vector relation k + p = q. These nonlinear interactions involve the 
neighboring Fourier components (k, p, q) that are excited in the local inertial range 
turbulence. The mode coupling interaction during the nonlinear stage of evolution 
leads to the formation of streamer-like structures in the magnetic field fluctuations 
associated with ky « 0, fc^; ^ 0. This is shown in Fig (1). Note that the streamer-like 
structures are similar to the zonal flows but contain a rapid variations, thus the 
corresponding frequency is relatively large. The temperature fluctuations in Fig (1), 
on the other hand, depict an admixture of isotropically localized turbulent eddies 
and a few stretched along the direction of the background inhomogeneity. We have 
performed a number of simulations to verify the consistency of our results in a 
strong turbulence regime. 

Figure (2) exhibits highly anisotropic nonlinear structures associated with mag- 
netic field vorticity (V^S) and generalized magnetic field vorticity {B — V^B). 
These structures, consistent with Fig (1), vary along the x direction and contain 
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Figure 3. Evolution of tlie total energy in tiie MEDV mode turbulence. A decay can 
be observed owing to the presence of small scale dissipation where much of the turbulent 
energy is concentrated. The inclusion of dissipation not only damps the smaller scales, but 
it also terminates the cascade processes that results in an enhanced inertial range MEDV 
mode spectrum. 

almost negligible gradients along the direction of inhomogeneity. The vorticities are 
increasingly dominated by nonlinearly generated fcj, « modes. The mode struc- 
tures in Fig (2) comprises elongated (along the y direction) smaller scale eddies 
because their spectra are propotional to and 1 + fc^, respectively. 



4. Evolution of energy in the MEDV mode turbulence 

The formation of relatively smaller scale fluctuations in the steady state MEDV 
mode turbulence (as displayed in Figs (1) and (2)) implies essentially that the 
inertial range turbulent spectrum is dominated by the higher Fourier modes. This is 
a characteristic of the forward cascade processes in that the energy containing large 
scales eddies transfer turbulent energy to the adjacent smaller scale until the process 
of energy migration is terminated by dissipative processes. To investigate the MEDV 
mode turbulent spectrum, we therefore include dissipation (arising from the plasma 
resistivity and electron thermal conduction involving electron-ion collisions) in the 
dynamical equations (5) and (6) to damp smaller scales and thus terminate the 
inertial range cascades. Note however that nonlinear interactions, in the absence of 
dissipation, conserves the total energy. 

It is worth emphasising here that the two-component (i.e. in B and T) system 
of equations describing 2D MEDV mode turbulence is characteristically different 
from the usual 2D drift wave or Euler turbulence systems, especially in terms of 
the nonlinear evolution. The presence of large scale structures have been routinely 
observed in the latter. One of the reasons is the inverse cascade processes that 
migrate spectral energy from the smaller to larger Fourier modes. By contrast, in 2D 
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Figure 4. Evolution of energy dissipation rates in the MEDV mode turbulence. During 
the inital decay, the dissipation rates vary rapidly followed by a steady state where decay 
rates are nearly negligible and follows a Kolmogorov-like energy cascade law. The latter 
is consistent with Fig (3) that depicts a constancy of energy in the steady state MEDV 
mode turbulence. 

MEDV mode turbulence the process of inverse cascade is inhibited conspicuously by 
the diagmanetic-like nonlinearity that depletes the large-scale flows. Furthermore, 
the interaction between the nonlinearity and the background inhomogeneous flows 
(associated dyT and dyB terms) leads to higly anisotropic nonlinear streamer-like 
strucutres (see Figs (1) and (2)). 

Since the energy is dominated by higher fc's in the MEDV mode turbulence, the 
inclusion of dissipation is expected to damp the turbulent energy. This is shown in 
Fig (3), where a significant decay of initial turbulent energy can be seen, which is 
followed by energy saturation in the steady state. The decay of turbulent energy 
is proportional to the dissipation. Thus larger is the dissipation, higher the energy 
decay rates are. After nonlinear interactions are saturated, the energy in the tur- 
bulence does not dissipate further and remains nearly unchanged throughout the 
simulations. Correspondingly, the energy transfer rate shows a significant change 
during the initial and nonlinear phases. However, when nonlinear interactions sat- 
urate, the nonlinear transfer of the energy in the spectral space amongst various 
turbulent modes becomes inefficient and the energy transfer per unit time tends to 
become negligibly small, as shown in Fig (4). 



5. Energy spectrum 

It is clear from the nonlinear structures associated with B, V^B and B — V^B that 
the mode structures are dominated by fcj, « streamer-like modes. Correspond- 
ingly, the inertial range spectrum should contain signature of the highly anisotropic 
streamer-like fiows. One thus expects that the resultant energy spectrum might de- 
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Figure 5. MEDV mode turbulence spectrum is much steeper than the usual inhomoge- 
neous drift wave turbulent spectrum in the forward cascade vorticity regime (fc-^). The 
observed spectrum in our simulation depicts a non Kolmogorov-like spectrum. The 
underlying steepness in the inertial range spectrum is mediated by the highly anisotropic 
nonlinear ly excited streamer-like structures for which ky ~ Q. 

viate from the standard Kolmogorov spectrum that is exhibited by the HM drift 
wave turbulence and its several variants. To investigate the spectral features of 
MEDV mode turbulence, we plot the inertial range energy spectrum in Fig (5). It 
is seen from this figure that the MEDV mode spectrum is much steeper than its 
usual 2D drift wave (isotropic) turbulence counterpart. The inertial range energy 
spectrum, in the latter, exhibits a Kolmogorov-like spectrum [18] in the vortic- 
ity cascade regime (see the adjacent curve in Fig (5)). The spectral index associated 
with the MEDV mode turbulence is close to fc~^ in the forward cascade vorticity 
regime. There is a greater disparity in the two spectra. While the disparity in the 
MEDV mode spectrum is caused essentially by the anisotropic turbulent cascades, 
it is interesting to know whether it can be understood from a Kolmogorov-like phe- 
nomenology. In view of pursuing this issue, we derive the MEDV mode turbulence 
spectrum based on the Kolmogorov-like dimensional arguments that are described 
in the following. 

The energy cascade is mediated essentially by the nonlinear convective forces in 
the MEDV mode turbulence and is proportional to kB^, where fc is a wavenumber 
associated with the magnetic field Bk in the spectral space. The corresponding 
energy transfer time can then be estimated as Tni ~ {kBk)~^ ■ The energy transfer 
rates, in the forward cascade vorticity regime, are determined by the vorticity per 
unit nonlinear time, i.e. 



e{k) 
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It should be borne in mind that the energy transfer rates, described above, are in- 
fluenced signiflcantly by the presence of the propagating MEDV modes. The linear 
normalized frequency of this mode can be determined from the dispersion relation. 
The linear interaction time of EDM with nonlinear eddies, in fc S> 1 dispersive 
limit, can be given as n ^ k/ky. This interaction time needs to be necessarily in- 
cluded in the energy transfer rates because the MEDV mode couples with nonlinear 
fluctuations. The modifled energy cascade rates are given in the following. 



'nl ky 

We next apply the Kolmogorov theory of turbulence that describes that the inertial 
range spectral cascade is local and depends on the wavenumber, as follows. 

Ek ~ e{k)°'k^. 

On substituting the energy transfer rates in the above expression, one obtains 

On equating the indices of identical bases, we get a = 2/3 and /? = —20/3. This 
further yields a Ek ^ fc"^*^/^ energy spectrum in the forward cascade vorticity 
regime. The spectrum obtained by including the MEDV mode interaction time in 
the Kolmogorov phenomenology is in close agreement with (if not exactly identical 
to) the observed energy spectrum (see Fig (5)). More steeper spectral index can be 
followed qualitatively from the entire energy spectrum that contains the anisotropic 
streamer-like flows. It can be described by the following expression. 

Ek ~ £(fc)2/3fc-20/3fc-2/3. (9) 

Thus, the deviation from the usual Kolmogorov-like k~^ spectrum increases with the 
anisotropy. The latter can be triggered either by the presence of the mean magnetic 
field or background inhomogeneous flows. The spectrum described by Eq. (9) depicts 
a much steeper spectral slope for streamer flows that have smaller magnitude of ky 
(or close to, but not entirely, zero). Note carefully that Eq. (9) does not reduce to 
a standard k~'^ spectrum describing the isotropic MEDV fluctuations, because the 
anisotropy of the spectral cascade is introduced in the energy cascade rates through 
a ratio of linear and nonlinear periods. This choice does not allow us to simply 
reduce the former to an isotropic turbulent spectrum. In this respect, the anisotropic 
turbulent spectrum described by Eq. (9) is only qualitative. A quantitative spectrum 
discerning the dependence on kx and ky may require an altogether different analytic 
approach and is beyond the scope of the present paper. 



6. Electron transport in magnetic structures 

We flnally estimate the electron transport in the presence of magnetic field struc- 
tures created due to nonlinear interactions between the MEDV modes. In the pres- 
ence of an ensemble of magnetic field structures the electrons perform a random 
walk and diffuse. Accordingly, there results and effective electron diffusion coeffi- 
cient, i'e/Zi which can be calculated from the cross correlation of the drift speed of 
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Figure 6. The MEDV mode turbulence transport is dominated by large scale streamer-like 
flows that are dominated by fcy ~ modes. The electron transport is higher initially and 
saturates in the steady state. 

the electrons in the turbulent magnetic field. We have 



where the angular bracket denotes ensemble spatial averages. The perpendicu- 
lar component of the electron fluid velocity associated with the MEDV mode is 
V_L = z X VB. Since the MEDV mode turbulence is confined in the 2D plane, the 
effective diffusion coefficient essentially describes the diffusion processes associated 
with the motion of electrons that produce fluctuating magnetic flelds in the x — y 
plane. We compute the evolution of the effective diffusion caused by nonthermal 
magnetic flelds of the MEDV modes. This is shown in Fig (6). The effective diffu- 
sion is larger during the early phase of the simulations. This is where the isotropic 
MEDV mode turbulence progressively evolve towards anisotropy and dominated 
increasingly by fcj^ « modes. Once the anisotropic streamers are formed, they 
begin to saturate such that no more energy transfer takes place. The corresponding 
steady state transport also shows saturation. This scenario is further consistent 
with the evolution of energy and the subsequent decay rates depicted in Figs (3) 
and (4), respectively. 

7. Discussion 

In this paper, we have used computer simulations to investigate the properties of 
2D MEDV modes in a nonuniform plasma containing the equilibrium density and 
electron temperature inhomogeneities. The equilibrium is maintained due to a bal- 
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ance between the electric force —neEox and the electron pressure gradient dpo/dx, 
where po = n{x)To{x). The electron temperature gradient can create a phase shift 
between the magnetic field and electron temperature perturbations, so that the 
MEDV mode grows provided that r] = dlnTo/dlnn > 2/3. Large amplitude MEDV 
modes interact among themselves and evolve in a turbulent state. The latter consists 
of magnetic filaments/streamers and localized electro heated regions in nonuniform 
plasmas. The resulting turbulent spectrum deviate from the usual Kolmogorov [18] 
or the Fyfe-Montgomery [19] scalings, which describe the complex fluid turbulence 
and drift wave turbulence in plasmas, respectively. Due to the random walk process, 
electrons would difi^use in the magnetic field structures that are created due to the 
nonlinear interactions between the MEDV modes. The magnetic field structures 
may also cause the electron heat transport on account of the electron temperature 
perturbation which are driven by the localized magnetic fields of the MEDV modes. 
The present investigation should help to understand the generation of flows, struc- 
tures and associated spectra producing anomalous transport of particles and heat 
in inertial confinement fusion plasmas. 
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